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Abstract The radionuchde ^^^Tb (Ti/2=6.89 d) is po- 
tentially important for internal radiotherapy. It is gen- 
erally produced through the is°Gd(n,7)isiGd ^^^Tb 
route at research reactors. In this work the possibility of 
its production at a cyclotron was investigated. Determi- 
nation of the excitation function of the ^^°Gd(d,x)^''^Tb 
production route and that of the disturbing ^^'^Gd(d,2n)- 
^^"Tb side reaction was done over the deuteron energy 
range up to 50 MeV using the stacked-foil technique 
and high-resolution 7-ray spectrometry. A comparison 
of this production route with the established (n,7) re- 
action at a nuclear reactor is made. 



1 Introduction 

Radioisotopes are essential for a variety of applications 
in medicine as diagnosis by scintigraphy and treatment 
of various diseases by internal radiotherapy. There is a 
rapid growth in the use of radionuclides for treatment 
of cancer in nuclear medicine paralleled by an increase 
of the diversity of the used radioisotopes. -'^^^Tb has 
appropriate decay properties for use in cancer therapy 
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— Its half-life of 6.9 days - long enough to allow trans- 
port to hospitals but short enough to avoid long 
term issues of waste handling after excretion from 
the patient. 

— Emits low-energy /3~- particles and low-energy Auger- 
electrons what results in a short cytotoxic range 
with minimal collateral damage to healthy tissue. 

— Emits a small amount of 7-radiation - enough to 
detect exactly where the radioisotope has been de- 
livered to. 

It has similar decay properties, biochemistry and in vivo 
property as the commercially available ^^^Lu radioiso- 
tope, but it emits more low-energy (Auger) electrons, 
meaning it could be better in treating smaller tumors. 
The radionuclide -'^^^Tb can be obtained in n.c.a. (non- 
carrier-added) form by irradiating highly enriched, mas- 
sive ^^'^Gd targets with intense neutron beams in nu- 
clear research reactors relying on the indirect ^^^Gd(n,7)- 
"'^^^Gd— T^^^-'^Tb production route [1]. This approach is 
similar to the indirect production route used for ^^^Lu. 
The thermal neutron cross-section of the ^™Gd(n,7) re- 
action (1.5 barn) is however lower than the one for the 
"6Yb(n,7)"6Yb^ I'^^fLu reaction (3 barn), therefore 
longer irradiation or higher flux is required for produc- 
tion of the same batch activity. The production possi- 
bility of using the photonuclear reaction ^®^Lu(7,p) on 
highly enriched -"^^^Lu was discussed in [S]. It will al- 
low to produce ^^-"^Tb with rather high specifc activity 
and low radionuclidic impurity, but the cross-sections 
for (7,p) reactions are known to be low. There is also a 
possibility to produce ^®^Tb through the {a,x) reactions 
from natural or enriched gadolinium [B] , but the excita- 
tion function reaches only 13 mb at 60 MeV. To increase 
the availability of ^^^Tb, we have investigated the pos- 
sible production by deuteron-induced reactions on Gd. 
No earlier investigations for this production route were 
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found in the literature. We should mention from the be- 
ginning that similar investigations were done for ^^^Lu 
by Hcrmanne et al. and by Manenti et al. [5] show- 
ing the advantage of the reactor production. We have 
used for our investigation gadolinium targets with nat- 
ural isotopic composition (^^^q^.q 2%, ^^^^0(1-2.18%, 
i55Gd-14.8%, i56Gd-20.47%, i57Gd-15.65%, i58Gd-24.84% 
i60Gd-21.86%) bombarded with a 50 MeV deuteron 
beam. As ^^^Tb can only be produced through reactions 
on the highest mass isotope ^^°Gd, among the numer- 
ous activation products assessed in this experiment we 
can single out cumulative production of ^^^Tb, via the 
direct (d,n) production and by the decay of ^^^Gd (3.66 
min, (3~: 100 % ). Results for the other investigated ra- 
dionuclides will be published separately. 

2 Experimental technique and data analysis 

The experimental techniques and data analysis were 
similar as described by us in the recent publications 
(c.f. [QlfTOllllj). Here we present shortly the most im- 
portant factors, and details, specific for this experi- 
ment. A first stack (series 1), containing "'^'Gd (83.9 
min), Sc target foils (105.27 min), aluminum degrader 
and monitor foils (27 /itm), was irradiated at the Cy- 
clone 90 cyclotron of the Universit Catholique in Lou- 
vain la Neuve (LLN) with a 50 MeV incident energy 
deuteron beam (42 min, 121 uA). A second stack (series 
2.) containing "'''Gd (83.9 m) target foils, aluminum 
degrader (98 fj,m) and Ti monitor foils (11 /im), was 
irradiated at the CGR 560 cyclotron of the Vrije Uni- 
versiteit Brussel (VUB) with a 21 MeV incident energy 
deuteron beam (60 min, 110 uA). The activity pro- 
duced in the targets and monitor foils was measured 
non-destructively (without chemical separation) using a 
high resolution HPGe 7-ray spectrometer. Three series 
of measurements at different detector-target distances 
were done starting about 2 h, 20 h and 200 h after EOB. 
The decay data were taken from the online database 
NuDat2 [T^. Effective beam intensities and the energy 
scale were determined by using the excitation func- 
tions of the 24Al(d,x)22'24Na and "''*Ti(d,x)48v reac- 
tions |13| simultaneously re-measured over the whole 
energy range. The figure of the re-measured excitation 
function in comparison with the IAEA recommended 
data [9lll3j can be found in our earlier publication on 
Sc, irradiated in the same experiment [11 . Concerning 
the two related isotopes (^^""^Tb and its ^^^Gd parent) 
we could not assess the short-lived ^^^Gd {T1/2 =3.66 
min), which has decayed out at the moment of the first 
series of measurements. Only the ^^^Tb 7-lines were 
identified. The 57.1 keV 7-line of ^^^Tb (see Table 1) 
was used for the cross-section calculations, as the 74.6 



Table 1 Decay characteristic of the investigated re- 
action products 



Nuclide 


Half- 


E-y (kcV) 




Contributing 


Q-valuo 




lifc 








(koV) 




6.89 d 


25.65135 


23.2 


-lB'JGd(d,n) 


458.3.96 


: 100 % 




48.91533 


17.0 










57.1917 


1.79 










74.56669 


10.2 








3.66 m 


102.315 


13.9 


-"S"Gd(d,p) 


3410.83 


/3- : 100 % 




314.92 


22.7 










360.94 


60.1 






i^lOTb 


72.3 d 


86.7877 


13.2 


-"^'■'Gd(d,2n) 


-3112.6 


fJ~ : 100 % 




197.0341 


5.18 










215.6452 


4.02 










298.5783 


26.1 










879.378 


30.1 










962.311 


9.81 










966.166 


25.1 










1177.954 


14.9 







7-line overlaps with the lead KX-rays and the 49 keV 7- 
line is contaminated with signal from the Tb, Gd and 
Eu KX-rays. At 25.6 keV 7-ray energy we had very 
low detector efficiency. Taking into account the low en- 
ergy of the used 7-line, the high Z and the thickness 
of the used target, correction was made for 7-ray self- 
absorption in the target (around 23%). For assessment 
of the contaminating long-lived ^^'^Tb many indepen- 
dent 7-lines are available (see Table 1). The uncertain- 
ties on the cross-section values are based on the stan- 
dard approach described in the ISO guide [HI and are 
obtained as the square route of the quadratic summa- 
tion of the independently contributing parameters. Un- 
certainty on the energy scale includes cumulative effects 
of initial beam energy, stopping estimation influenced 
by target thickness and beam straggling. 

3 Theoretical calculations 

The cross-sections of the investigated reactions were 
calculated using the pre-compound model codes ALICE- 
IPPE [H] and EMPIRE-II [W . To improve the descrip- 
tion of the experimental data of deuteron induced reac- 
tions at lower energies modified model codes were used. 
In the ALICE IPPE-D and EMPIRE-D code versions 
the direct (d,p) channel is increased strongly, which re- 
flects in the results for all other reactions in both codes 
[ITlfTS] . A comparison with the data retrieved from the 
online TENDL 2011 and TENDL 2012 libraries, based 
on the latest version of the TALYS code system [19] is 
also made. We presented both TENDL 2011 and 2012 
versions in the figures to see the improvement of the 
description. 

4 Results 

When supposing that highly enriched "'^^"Gd targets are 
used, the measured production cross-sections on "°*Gd 
can be transformed into reaction cross-sections on pure 
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Table 2 Experimental 
i6"Gd(d,x)ieiTb(cum) and 
tions 



cross-sections for 
i60Gd(d,x)i60Tb reac- 



I '«°Gd(d,x)«'Tb(cum) 



E ±ziE (McV) 


Cross-scction(cr)±zicr (mb) 


1 


^•^-^TbCcum) 1 ^*^*'Tb 


Series 1 


49.2 


0.3 






49.6 


8.4 


47. S 


0.3 






52.9 


13.0 


46.4 


0.3 


76.1 


15.9 


52.2 


6.3 


45.0 


0.4 


66.1 


36.8 


57.5 


14.8 


43.6 


0.4 






52.2 


10.1 


42.1 


0.4 






45.4 


11.5 


40.5 


0.4 


77.0 


8.3 


48.7 


12.4 


39.0 


0.5 






65.8 


14.2 


37.4 


0.5 


72.2 


29.0 


72.4 


13.6 


35.7 


0.5 


68.4 


7.4 


74.2 


19.1 


34.0 


0.5 






100.0 


17.8 


32.2 


0.6 






79.6 


15.0 


30.4 


0.6 


105.2 


11.3 


109.4 


17.3 


28.4 


0.6 


104.2 


11.2 


107.8 


16.8 


26.4 


0-7 


112-5 


39.5 


105.5 


16.2 


24.2 


0-7 


147-0 


15.9 


132.4 


18.4 


21.9 


0.7 


159.9 


17.3 


147.1 


19.8 


19.4 


0.8 


102.5 


20.9 


197.2 


23.3 


16.6 


0.8 


203.7 


22.1 


329.4 


37.4 


13.5 


0.9 


223.1 


24.2 


553.1 


60.8 


Series 2 


19.8 


0.3 


144.8 


26.0 


180.3 


20.7 


18.3 


0.3 


165.6 


22.3 


239.0 


26.9 


16.6 


0.4 


176.4 


22.9 


334.6 


37.2 


14.9 


0.4 


180.3 


21.0 


432.1 


47.1 


13.0 


0.4 


249.7 


28.8 


663.7 


72.3 


10.9 


0.5 


292.9 


33.4 


544.0 


59.9 


8.4 


0.5 


174.0 


19.6 


378.9 


41.6 


5.4 


0.6 


5.1 


0.9 


14.1 


1.8 



^^"^Gd. From the point of view of radionuclidic purity 
only the long-hved ^^"Tb (Ti/2=72.3 d) should be con- 
sidered as a possible and relevant contaminant up to 
the threshold of the is°Gd(d,6n)i56Tb reaction ( Q = 
-33560.68 keV). Other Tb radioisotopes produced at 
lower energy are either stable (-'^^^Tb) or have very long 
half-life (158Tb (150 a) and ^^Tb (150 a)). The cumu- 
lative excitation function for the i^*'Gd(d,x)i^iTb reac- 
tion is the sum of i60Gd(d,n)isiTb and i^"Gd(d,p)i6iGd-^ 
i*5iTb and the excitation function of the if^°Gd(d,2n)if^°Tb 
reaction, compared to the results of three theoretical 
codes are presented in Fig. 1 and Fig. 3. The numerical 
data are presented in Table 2. 

The uncertainties on ^^^Tb cross-sections are large 
(and the data are scattered), due to the low statistics 
of the 57.2 keV 7-line and to the uncertainty of the 
peak fitting in a complex spectrum. The difference be- 
tween the particular uncertainties can be explained by 
the fact that different spectra were measured after dif- 
ferent cooling times and for different measuring times, 
so some of them have worse statistics than others. 

In Fig. 1 we have included the results for cumula- 
tive production of ^^^^Tb using the phenomenological 
systematics of (d,p) reaction taken from contribution 
of A.V. Ignatyuk to FENDL-3 Activation Library [20] 
(Fig. 1, sys) The most striking information that can be 
derived from the Fig. 1 is: 

~ The direct contribution to ^^^Tb production is very 
low. It is based on the well-known systematics of 
the (d,n) reactions in this mass region and from our 
theoretical calculations. 

— The Tb is produced mostly indirectly. 




o this work (cum), s 
• this work (cum 

ALICE-D (cum 

EMPIRE-D(ci 

— TENDL 2011 (CL 
■ - TENDL 2012 (cu 

Sys (oj 

■■■■ ALICE-D (dir) 
EMPIRE- D ( dir) 
TENDL 2011 (dir) 
TENDL 2012 (dir) 



Deuteron energy (MeV) 

Fig. 1 Excitation function for the reactions 
^®'^Gd(d5x)^^^Tb(cum) and comparison with the theo- 
retical calculation 



^"Gcl(d,x)i6iTb 




o 161Tb exp, this work, series 
• 161Tb exp, tiiis work, series 
a177gLuexp Manenti 
■ 177gLu exp. Hermanne 
+ 1 77gLu exp. Tarkanyi Unp. 



Deuteron energy (MeV) 

Fig. 2 Comparison of tlie excitation 
for the reactions isoq^jj-jJ x)ieiTb(cum) 
i^6Yb(d,x)i'^'^!'Lu(cum) 



functions 
and and 



~ The agreement with the results of the modified AL- 
ICE and EMPIRE codes is acceptable good. 

- Both TENDL 2011 and TENDL 2012 underestimate 
significantly the cumulative cross-section (and con- 
sequently the (d,p) by a factor of 2). The agreement 
with the TENDL 2012 resuhs is shghtly better. 

The cumulative excitation functions of the ^^°Gd(d,x) 
^^^Th and the "6Yb(d,x)i^7sLu [HEIIIII are similar 
both in shape and in magnitude as it is illustrated 
in Fig. 2. (The direct and indirect production of the 
i77m]^^j is negligible due to the high spin value J'^ = 
23/2^). In Fig. 2 our new data (preliminary, unpub- 
lished) measured for ^^^Yb(d,x)^^^f Lu reaction up to 
40 MeV are also shown for comparison at the high en- 
ergy range. The cumulative excitation function for the 
^^"Gd(d,2n)^^°Tb reaction is shown in Fig. 3 in com- 
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Deuteron energy (MeV) 



Fig. 3 Excitation function for tfie reactions 
^®°Gd(d,2n)^^''Tb and comparison with tfie theoretical 
results 



parison with the results of the model codes and the 
following conclusions can be drawn: 

— The contaminating ^^'^Tb is produced directly via 
i*^°Gd(d,2n) reaction. 

— The agreement with the results of the modified ALICE- 
D and TENDL 2011 and TENDL 2012 is acceptable. 

— The results of EMPIRE-D overestimate the experi- 
mental values. 

— At high energies above 20 MeV the results of all 
model codes are systematically too low. 

— The maximum of the (d,2n) cross-section in this 
mass region is high, in agreement with the system- 
atics. 



Deuteron energy (MeV) 

Fig. 4 Integral yields for i80Gd(d,x)i''iTb and 
i'30Gd(d,x)i''0Tb reactions 



Table 3 Comparison of production routes of ^^'^Th 





Reactor 


Accelerator 


Reaction 


J-"^^Gd(n,7)^*^J-Gd^-L'^J-Tb 


J-*^'JGd(d.x)-L*^^Tb 


Flux (part/sec 


1.5*10-'-^ 




1.53*10-'-^p/scc 
(250 mA) 


6.1*10-'-^p/sec 
(1000 mA) 


Energy range 


thermal 


thermal 


30-5 MeV 


30-5 MeV 


Target mass (mg) 


40 


40 


3 


3 


Average cross- 
section (mb) 


1500 


1500 


125 


125 


Yield ratio to high 
flux reactor 


1 


1/50 


1/266 


1/66.5 


Irradiation 


multichannel 


single channel 


Chemical separa- 
tion 


yes 


yes 


Enrichment 


highly enriclicd 


liighly enriched 


Specific activity 


higher 




Disturbing reaction 






-^'^"Gd(d,2n)^^l'Tb 



also at accelerators. The main parameters for compar- 
ison of the reactor and accelerator routes are collected 
in Table 3. 



5 Thick target yields 

From fitted curves to our experimental cross-sections 
integral thick target yields were calculated for produc- 
tion of ^^^Tb(cum) and ^^"Tb. The results for physical 
yields [12] are presented in Fig. 4. No earlier experimen- 
tal thick target yield data were found in the literature 
for comparison. 

6 Comparison of the accelerator and reactor 
production routes 

The presently used route for production of carrier free 
i6iTb is the indirect i60Gd(n,?)i6iGd-^ i^^Tb process 
at high flux nuclear research reactors or at spallation 
neutron sources [TJ[23] followed by a chemical separation 
of radionuclide of interest from the target material. Ac- 
cording to the present study the ^^^Tb can be produced 



6.1 Production yield 

The Table 3 shows that from the point of view of pro- 
duction yield of ^^-'^Tb the recent accelerators are not 
competitive with a high flux reactor. The productivity 
of a high intensity accelerator is close to the productiv- 
ity of standard reactors. 

6.2 Targetry 

Use of charged particle beam requires special, high power 
targetry, which is not a simple question. The thermal 
conductivity of Tb is 40 times lower than that of Cu 
(Tb-10.6 Wm^Ki and Cu-400 Wm^K^). Using highly 
enriched target material, different target preparation 
methods can be considered. Electrodeposited in metal 
form can be used taking into account that the Gd metal 
is relatively stable in dry air [53]. Different compounds 
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Radionuclide impurity 




- 0.1 

— I hour (number) 
• 24 hour (number) 
-24hDur(acrivitv) Q QS 

- 1 hour (activity) 



Deuteron energy (iVIeV} 



Fig. 5 5. Radionuclidic impurity of the irradiated highly en- 
riched ^®'^Gd after different irradiation times. Both the ratio 
of the number of radionuclides (left) and the ratio of activities 
(right) are presented. 



(oxide, chloride) can be used, prepared by pressing or 
in molten form. 



6.3 Radionuclidic purity, specific activity 

In case of neutrons, in principle there is no disturbing 
impurity by using highly enriched targets followed by 
an effective chemical separation. By using deuterons up 
to 40 MeV of commercial cyclotrons on "°*Gd targets a 
considerable amount of I53,i54,i55,i56,i57,i58,i60rpi3 

be produced. Therefore the radionuclidic impurity will 
be high. By using highly enriched -'^^"Gd (100%) up to 
20 MeV bombarding energy, only ^^"Tb and ^^^Tb ra- 
dioisotopes will be produced. The energy domain of the 
excitation functions of the ^^^Tb and -"^^^Tb is nearly 
similar but the cross-section is 2-3 times higher for the 
production of ^^°Tb. Therefore it is difficult to reduce 
the radionuclidic impurities and to increase the specific 
activity by using an adapted energy domain (see Fig. 5). 
Roughly two times more ^^°Tb than ^^^Tb nuclei will 
be present in the irradiated sample, even using short 
irradiations (see Fig. 5). The ratio will even be worse 
for longer irradiations (more than 100 h), because it 
will increase with decay after EOB approximately lin- 
early. Because the allowed radionuclidic impurity may 
not exceed 1-2% in the case of human application, this 
way of production is completely unsuitable for medical 
purposes. 

The conclusion based on the short discussion above 
is that accelerator production of ^^^Tb is not compet- 
itive with the reactor production both in the aspect 
of commercially attractiveness (investment in irradia- 
tion facility not considered) and product quality. The 



only advantage of ion-accelerators production could be 
the easier availability in small batches for research pur- 
poses. 



7 Comparison of the experimental and 
theoretical cross-section data 

To model a reliable and accurate theoretical description 
of deuteron induced activation mechanisms knowledge 
of experimental data is very important. The interaction 
of the deuteron with the target nuclei is very complex, 
as it includes the deuteron break up process and the 
related stripping and pick-up mechanisms, statistical 
emission, etc. Possibilities to take these different pro- 
cesses into account in model codes have been proposed 
in different-simple and complex ways [TSJUSKIB] . The 
comparison with the present experimental data also 
shows that the theoretical description still needs signif- 
icant improvement in the TALYS code, especially for 
(d,pxn) processes. In case of the modified ALICE and 
EMPIRE codes the agreement is better, but it is a result 
of a correction based on a phenomenological simulation 
of direct (d,p) and (d,t) transitions, including a "re- 
duction factor" of the compound nucleus cross-section 
due to direct processes based on systematic behavior of 
the experimental data. Our systematical study of the 
deuteron induced reactions [27,28 can contribute effec- 
tively to further improvement of the theoretical treat- 
ment. 
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